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Introduction 


A two-year investigation of the trailing vortex shed by a rectangular 
NACA 0012 wing is in progress under the above grant. Tlie objectives of 
this investigation are to observe the turbulence structure and spectral 
characteristics of this flow over a range of conditions and to incorporate 
these observations into the blade-wake interaction (BWI) noise-prediction 
method of Glegg(1989). This report is divided into four sections as 
follows :- 

1. Measurements performed during the first year. 

2. Presentation and discussion of a representative sample of the results. 

3. Implications for the BWI noise prediction method. 

4. Re-evaluation of work planned for the second year. 

1. Measurements performed during the first year. 

Experiments are being performed in the Virginia Tech Stability 
Wind Tunnel. So far the measurements have included helium-bubble flow 
visualizations and single-point hot-wire velocity measurements. 

The visualizations, which are described in detail by Devenport and 
Sharma (1990), were taken over the range of Reynolds numbers and angles 
of attack listed in table 1. These showed the vortex to be relatively stable 
(moving less than ~5% chord across the tunnel) and insensitive to probe 
interference, making useful hot-wire velocity measurements possible. 

The following velocity measurements have been made at conditions 
listed in table 2. 

(a) Single hot-wire measurements at the wing trailing edge. These 
include detailed profiles and power spectra of the streamwise 
component of velocity. The purpose of these is to reveal the 
properties of the boundary layer shed from the wing. 

(b) X-array hot-wire measurements in the vortex. Detailed profiles of 
all mean velocity and Reynolds stress components. Single-point 
spectra of all velocity components at representative radial positions. 

(c) Triple hot-wire measurements in the vortex. Detailed profiles of all 
mean velocity and Reynolds stress components. Single-point spectra 
of all velocity components at representative radial positions. Long- 
time-record measurements at selected radial locations to deduce the 
statistical characteristics of the vortex motions. 



Z Presentation a nd Discussion of Experimental Results. 

Figure 1 shows the coordinate directions (X,Y,Z) (U,V,W) to be 
used in presenting results. The wing chord c, of 8", will be used to 
normalize distances. The nominal undisturbed free-stream velocity U re f 
will be used to non-dimensionalize velocities. The thickness and state of the 
boundary-layer leaving the wing are listed in table 3 for each of the flow 
conditions. 

Figures 2, 3 and 4 show velocity measurements made during a Z- 
wise traverse through the center of the vortex at X/c = 30, Re c = U re fc/v = 
400000, with the wing at 5° angle of attack (rotating from the X to the Z 
axis). Note that Z is measured relative to the center of the vortex. 

In the vicinity of the vortex center (figure 2) the mean-velocity field 
is much as would be expected given the results of Mason and Marchman 
(1972). The tangential velocities associated with the vortex are clearly 
visible in the W component profile (figure 2(b)). These increase to a peak 
at the edge of the vortex core, which appears to have a diameter of about 
0.72",(=0.09c). The small the core size is circumstantial evidence for the 
relative stability of the vortex and its insensitivity to probe interference. 
The U profile (figure 2(a)) shows an axial velocity deficit of about 
0.12U re f. The normal stress profiles (figure 2(c)) are dominated in this 
region by a strong central peak. Here v 2 /U re f 2 and w 2 /U re f 2 reach peak 
values of 0.0064 (8% turbulence intensity) and u 2 /U re f 2 a peak of 0.0007 
(2.6%). The large values of v 2 and w 2 are almost certainly due to a 
combination of the the small lateral motions of the vortex observed in the 
flow visualizations and the steep W gradient in the core. Approximate 
modeling of these motions, described in the appendix, show them to have 
an rms amplitude of less than 2mm. Observations of the velocity signals 
during the measurement suggest they have a frequency between about 1 and 
10 Hz. 

Away from the vortex center (figure 3) the results reveal a flow 
structure somewhat different from what we had expected. Here the 
measurements show two, if not three, distinct half turns of the wing wake 
as it spirals around the vortex. These appear as peaks in the normal and 
shear-stress profiles (figures 3(c) and (d)) and inflections and depressions 
the W and U profiles (figures 3(b) and (a)) respectively, centered at Z = 
+6", -2.5" and +1.5". In contrast to the results and discussion of previous 
workers (in particular Phillips and Graham (1984)) there appears to be 
little or no region surrounding the core in which the successive turns of the 
wing wake have merged to form a continuous axisymmetric structure. The 
behavior of the turbulence stresses in the spiral wake is especially 


interesting, (the upwash component) is easily the largest stress, being 
typically 1.5 times u^ and twice v^. The magnitude of the normal stresses 
fall by about 40% on each half turn of the spiral. The shear stresses (figure 
3(d)), of which uw is the largest drop, even more rapidly. This suggests 
that the greater circumferential shear experienced by the wing wake 
towards the vortex center inhibits the development or motions of large 
stress -producing turbulent eddies. In interpreting the finer details of the 
turbulence structure it is important to realize that there is no significant 
mean motion along the spiral towards the center of the vortex and that, 
after a few chord lengths of initial development just downstream of the 
wing, the tangential velocity field probably changes little in scale or 
structure with stream wise distance. In other words turbulence structure 
appearing at a particular radial position at X/c=30 presumably developed at 
that position under the action of an fairly constant mean velocity field. 

Figure 4 shows autospectra of upwash velocity (w) measured 
different distances from the vortex center (at Z = -4”, -2", -1", -0.7", - 
0.35" and 0") normalized on w^. Note that the discrete peaks in these 
spectra at frequencies above 2kHz are a result of electrical noise. At Z = - 
4" and -2" (to either side of the 2nd half turn of the spiral wake, figure 
3(c)) the spectra are almost identical despite the fact that the turbulence 
levels at these two points are quite different. These spectra have a broad 
peak centered at about 300Hz, corresponding to a streamwise lengthscale 
(assuming Taylor's hypothesis) of about 4", a distance approximately equal 
to the width of the wing wake at X/c = 30. The existence of such large 
structures seems improbable, however, when one considers the significant 
tangential-velocity gradient across the wake at this location (figure 3(b)). 
Moving in towards the vortex center the upwash spectrum changes rapidly. 
Most noticeable is a large increase in the energy at very low frequencies 
(presumably associated with the lateral motions of the vortex) and the 
development of a strong peak around 600Hz, yet to be explained. 

Because of the possible importance of the spiral structure to BWI 
noise generation we have included normal stress profiles showing its 
behavior over a range of conditions in figure 5. The spiral widens with 
distance downstream of the wing and the turbulence levels within it fall 
(figures 5(a), (b) and (c)). As one would expect, increasing the angle of 
attack, and thus the vortex circulation (figures 5(c), (d) and (e)) has exactly 
the same qualitative qualitative effect as moving downstream. Changing the 
Reynolds number (figures 5(c), (f), (g) and (h)) has some effect on 
turbulence levels within the spiral wake but does not otherwise appear to 
influence the flow structure. 



3. Implications for the B WI noise prediction method 

The BWI noise prediction method was based on the measurements of 
Phillips and Graham(1984) which showed the turbulence intensity decaying 
monotonically from the center of the vortex. A simplified model was 
developed for the distribution of the turbulence and this gave a reasonable 
fit to the measured noise spectra, but spectral shape was incorrect. In this 
section the original BWI turbulent flow model will be compared with the 
current set of measurements and suggestions made as to how this model 
could be improved. 

For the flow conditions at a Reynolds number of 400000 the original 
BWI model predicts the distribution of normal turbulent stress shown in 
figure 6. This identifies some important departures from the measured 
levels given in figures 2(c) and 3(c). First note that the peak turbulent 
stress is predicted as 2.6x1 0 4 in comparison with the measured peak level 
of 64xl0 -4 . The large measured levels in the vortex core have been 
attributed to the lateral motion of the vortex whereas the BWI prediction 
scheme assumes purely turbulent motion relative to a fixed vortex. The 
significance of the large measured velocity fluctuations is questionable 
since it is argued that the lateral motions of the vortex do not cause BWI 
noise for the following reasons. First the velocity fluctuations only occur 
close to the core which is unlikely to be split by the blade. Secondly the 
vortex core is symmetric and so there will be equal upwash and downwash 
on the blade at any instant, causing an inefficient mechanism for producing 
an unsteady load. Thirdly the frequency spectra of the lateral vortex 
motions show more energy at low frequencies than would be expected 
from BWI noise measurements. We conclude therefore that, although the 
measured turbulence levels are highest close to the center of the vortex due 
to lateral motions, these are not significant for BWI noise. 

Better insight is obtained by comparing figure 6 with figure 3(c) 
which is an expanded plot of figure 2(c). The measured and predicted 
turbulence levels, outside of the central core, are of the same magnitude 
although their distributions are quite different. It is clear that the original 
BWI prediction method uses an inadequate model of the turbulence 
structure and further consideration must be given to modelling this more 
accurately. 

In the previous section the spiral structure of the flow about the 
vortex core was discussed. If we assume a flow with vorticity concentrated 
in the vortex core, which is not completely accurate for this flow, then a 
simple model of the spiral can be developed which enables some 
conclusions to be drawn about BWI noise generation. The azimuthal 
velocity outside the vortex core can be described as 
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where y is an empirical constant, chosen as 7^0.8, and 0 is the momentum 
thickness of the wake. The upwash spectrum is then calculated as 


G ww( f ) = 


w 2 L 

tcU 


l + (14/3)kj 

l(l + k?) U/6 ] 


where 

- coL 

k, = 

1 1.34U 

Spectra based on this model are shown in figure 8 and clearly indicate 
more low frequency energy and a slower decay at high frequencies than 
measured in the wake spiral. These features were also found in the 
predicted BWI noise spectrum when they were compared with measured 
noise levels. This suggests that if the turbulence spectra measured here 
were used for noise predictions better results would be obtained. 

4. Re-evaluation of Work Planned for the Second Year. 

The BWI noise prediction model needs to be updated to incorporate a 
suitable model of the flow in the wake spiral. This will require: 

-models for the spiral wake boundaries, and the turbulence stress 
distributions, 

-models of the spectral distribution and lengthscales of the turbulence. 
To develop these models further measurements are clearly needed. 
Specifically, 

-detailed measurements of turbulence stress distributions throughout 
the spiral, 

-measurements showing the variation of velocity spectra along and 
across the spiral, 

-measurements of lengthscales along and across the spiral at a 
number of different positions within the spiral. 

In addition to providing quantitive input for the BWI noise 
prediction method, these data are needed to obtain a clear physical 
understanding of the flow. Such an understanding is essential if the models 
to be developed for the noise prediction method are to have any generality. 
With this in mind we beleive it will be most effective to concentrate on 
obtaining a detailed and complete data set at one flow condition (say 
X/c=30, Rec=400000, 5° angle of attack). It should be possible to infer the 
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The Unsteady Velocity Components 


Consider the vortex center as displaced from the center of the co- 
ordinate system by z which is a random function of observer time at a 
fixed streamwise location (see figure 3). A turbulent eddy located at x 
relative to the origin of the co-ordinate system, is displaced from the center 
of the vortex by y, such that x=y+z(t). The mean velocity field of the 
vortex is defined as V(y) and the unsteady turbulence component is v(y,t) 
and are a function of the displacement from the vortex center. The velocity 
measured by a fixed probe will have a mean component U(x) and an 
unsteady component u(x,t), given by 

U (x)+u(x,t)= V (x-z(t))+v(x-z(t),t) 

If the displacements of the vortex are small in comparison with the the core 
radius then we can expand the right hand side of this equation in a Taylor 
series, to first order, 


U(x) + u(x,t) = V(x) + v(x,t) - (z(t).V)(V + v) 

This shows that the measured mean component U(x) is, to first order, the 
same as the mean flow about the vortex. The unsteady component is a 
combination of the turbulence associated with the vortex and velocity 
fluctuations induced by the time varying lateral displacements. 

We assume that v and z=(0,oc,p) are uncorrelated so that the mean 
square velocity components are 



If the vortex displacements can be assumed to be isotropic with an rms 
displacement e, so that 


<cc 2 >=<p 2 >=£ 2 and <ap>=0 (A4) 

the three components of the normal turbulent stresses due to the lateral 
vortex motion can be defined as 




where V=(U,V,W). By using the mean velocity profiles described in the 
previous section we obtain 
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The measurements were taken on y=0 and so these equations can be 
simplified to give the three components of the normal turbulent stresses as 




By using equations A1 and A2 we obtain 
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Notice here that as 'r' becomes small we can approximate 
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In figure 2(c) we see that <v’ 2 >=<w ,2 >=0.0064U re f 2 and since v o =0.24U re f 
and r o =0.36'\ we can compute e as 0.0683"(1.73mm). 

Figure A4 shows the predicted normal turbulent stresses obtained 
using this model and the agreement with the measurements in Figure 2(c) is 
excellent. The only noticeable error is in the predicted value of <u' 2 > close 

to z -0. However the effects of finite probe size and low turbulence levels at 
this point may cause he predicted dip to be filled in. 

The turbulent shear stresses can also be predicted using this 
approach, since 6 


<uv)= B 4feY^U^Ysvi 

Iv AdyJ v dz A dz j 


(u'w') = e : 


1 + f— Y— ll 

V dz A dz 


v a y A 9 y J 


Inspection of equations A5,A6, and A7 shows that the only non-zero 
turbulent shear stress on y=0 is the <u'v’> component which is given by 
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In terms of the simplified model this can be defined as 
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and is plotted in figure A5. Comparing figure A5 with figure 2(d) shows 
that there is a significant discrepancy between the predicted and measured 
levels. First the predicted values if <uv> are much higher than measured 
and secondly the measured values of the <wv> shear stress is not 
negligible. These features of the flow do not appear to be explained by this 
simplified model 


Reynolds number 


Angle of attack 


1.3 x 105 2.5° 

5° 

7.5° 

2.6 x 105 5° 

7.5° 

4.0 x 105 2.5° 

5° 

7.5° 

5.3 x 105 5° 


Table 1. Flow conditions for helium bubble visualizations (120 grade 
sandpaper trip). 


Reynolds no. 

Angle of attack 

Axial pos 
z/c 

Boundary layer 
trip 

1.3 x 105 

5° 

30 

Glass beads 

2.6 x 10 3 

5° 

30 

Glass beads 

4.0 x 10 3 

0° 

30 

Glass beads 


2.5° 

30 

Glass beads 


5° 

30 

Glass beads 


5° 

25 

Glass beads 


5° 

20 

Glass beads 

5.3 X 105 

7.5° 

30 

Glass beads 

5° 

30 

Glass beads 


5° 

30° 

None 


Table 2. Flow conditions and locations of hot-wire measurements. 



Angle of attack 

Side 

Reynolds number 

Thickness 1 

State 

(degrees) 


Re c 

(inches) 


0 

S 

400000 

0.353 

Turb. 


P 


0.357 

Turb. 

2.5 

s 

400000 

0.390 

Turb. 


p 


0.323 

Turb. 

5.0 

s 

130000 

0.420 

Trans. 


p 


0.270 

Trans. 


s 

260000 

0.413 

Turb. 


p 


0.300 

Low Re Turb. 


s 

400000 

0.418 

Turb. 


p 


0.292 

Turb. 


s 

530000 

0.437 

Turb. 


p 


0.298 

Turb. 

7.5 

s 

' 400000 

0.552 

Turb. 


p 


0.270 

Turb 

Table 3. Boundary-layer thickness and state at the 

wing trailing edge 

with the glass-beads 

trip. - suction side, (P) - 

pressure side. State infered from mean-velocitv and turbulence 

intensity profiles. 






‘Defining the boundary layer edge as the point of 2% turbulence intensity. 
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Shear Turbulence Stresses Normal Turbulence Stresses 



(c) Profiles of the normal turbulence stresses. 



(d) Profiles of the turbulence shear stresses. 

Figure 2. Velocity measurements near the vortex center at X/c = 30, Re c = 400000, 5 
angle of attack. 




V/Uref, W/Uref 




Figure 3. Velocity measurements away from the vortex center at X/c - 30, 
Re c = 400000, 5° angle of attack. 






Figure 3. Velocity measurements away from the vortex center at X/ c 
Re c = 400000, 5° angle of attack. 
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Normal Turbulence Stresses Normal Turbulence Stresses 




Normal Turbulence Stresses Normal Turbulence Stresses 




Figure 5. Normal turbulence stress profiles measured away from the vortex center. 



Normal Turbulence Stresses Normal Turbulence Stresses 



(e) X/c = 30, Re c = 400000, 2.5° angle of attack. 



(f) X/c = 30, Re c = 130000, 5° angle of attack. 


Figure 5. 


Normal turbulence stress profiles measured away from the vortex center. 
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Figure 5. Normal 
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Figure Al: The axial velocity profile predicted using equation A2 with 
UD=0.12Uref and r o =0.36 M 
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Figure A2: The azimuthal velocity profile predicted using equation 
Al with v o =0.24U re f and r o =0.36 M 




Figure A3: The co-ordinate system relative to the center of the vortex. 



motion model. 



Turbulent Shear stress <uv>/Uref A 2 



Figure A5; The predict 



